Microscopic theory of magnetism in Sr3lr207 



Jean-Michel Carter^ and Hae- Young Kee^'^'j^] 

^Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7 Canada 
^Canadian Institute for Advanced Research, Toronto, Ontario Canada 

An intriguing idea of spin-orbit Mott insulator has been proposed to explain magnetic insulating 
behavior in various iridates. This scenario relies on the strength of the spin-orbit coupling being 
comparable to electronic correlations, and it is not a priori obvious whether this picture is valid for 
all iridates. In particular, Sr3lr207 exhibits a small charge gap and magnetic moment compared 
to Sr2lr04, questioning the validity of such hypothesis. To understand the microscopic mechanism 
for magnetism in Sr3lr207, we construct a tight binding model taking into account the full t2g- 
orbitals, the staggered rotation of the local octahedra, and the bilayer structure. The bands near 
the Fermi level are mainly characterized by the total angular momentum Jeff = 1/2, except below 
the r point, supporting a reasonably strong spin-orbit coupling picture. A first order transition to 
a colhnear antiferromagnet via multi-orbital Hubbard interactions is found within the mean field 
approximation. The magnetic moment jump at the transition is consistently smaller than Sr2lr04, 
originated from the underlying band structure of a barely band insulator. Given the small charge 
gap and moment observed in Sr3lr207, the system is close to a magnetic transition. A comparison 
to a spin-model is presented and connection to the Mott insulator is also discussed. 



I. INTRODUCTION 

The concept of Mott insulators is prevalent in strongly 
correlated materials, referring to a strong electronic cor- 
relation driven insulator that violates the conventional 
band theory. Often, Mott insulators are accompanied 
by antiferromagnetic ordering, but genuine Mott insula- 
tors have a robust charge gap above the Neel temperature 
where the magnetic ordering disappears. Mott insulators 
have been found in correlated electronic systems with 3d- 
orbital materials including the high transition tempera- 
ture cuprates, and are considered seeds of exotic phases 
when doped by holes or electrons. However, with heav- 
ier atom systems, the outer shell electron wavefunctions 
are less localized leading to weaker electronic correlation, 
and thus band theory should be a good starting point to 
model the behavior of electrons in such solids. 

Surprisingly, Sr2lr04, a material with heavy 5d Ir 
atoms, exhibits a magnetic insulating state, despite hav- 
ing partially filled 5d-orbitals. One missing ingredient 
that is relevant in these heavy elements is spin-orbit 
coupling (SOC). In particular, when the SOC strength 
is comparable to that of electronic interactions, un- 
derstanding their interplay becomes challenging. Iri- 
dates with 5d-orbitals offer such a playground to inves- 
tigate their combined effects. It has been suggested that 
due to the narrowing of the bandwidth induced by the 
strong SOC, the effect of Hubbard interactions is am- 
plified, leading to an insulating state in some layered 
perovskite^ ^ and pyrochlore^° ir idate s. This state was 
called a spin-orbit Mott insulator.'^'^ While this pro- 
posal of Mott i nsulat or for Sr2lr04 itself was challenged 
in recent studieJ^^EH^ the idea ha s been quickly applied to 
other iridates including Sr3b:20'} ^^ ^ ^^ \ its bilayer sister. 

However, optical conductivitj^l and transportPE^ data 
have shown that the charge gap in the bilayer system is 
significantly smaller than in the single layer compound. 
In addition, the magnetic structures in these two sys- 



tems differ; Sr3lr207 displaying a collinear antiferromag- 
netic structure with moments aligned with the c-asdsr^^^ 
whereas Sr2lr04 is well-known to have a canted antifer- 
romagnetic structure with moments in the ab-plan d^ ^ ^ l ^ l 
This may be related to the crystal structure of the 
two compounds. Sr3lr207 crystalhses in a Bbcb space 
groupP^l whereas Sr2lr04 does with an I4i/acd space 
groupP. Another major difference between these two 
materials is the giant magnon gap seen in Sr3lr207fSl 
whereas in Sr2lr04, the magnon gap, if it exists, is too 
small for detection with current resonant inelastic x-ray 
scattering (RIXS) resolutiorF^. Based on these differ- 
ences, it is not clear whether a strong SOC approach is 
valid Sr3lr207. 

In this paper, we build a tight-binding Hamiltonian 
for the bilayer Sr3lr207 taking into account the full 
t2g manifold, the local staggered rotation of the octa- 
hedra, and bilayer coupling, to understand the validity 
of strong SOC picture and the interplay among SOC, 
electronic correlation, and the crystal structure. The pa- 
per is organized as follows. In Sec. |llj using a Slater- 
Koster theorjl^l^ the tight binding parameters are es- 
timated. Comparing this band structure with r ecent 
angle-resolved photoemission spectroscopy (ARPES)^^^ 
measurements and first-principle calculation^, we deter- 
mined the SOC strength. Indeed the bands near the 
Fermi level are mainly Jeff = 1/2, except near the T 
point below the Fermi level, favoring a reasonably strong 
SOC limit. Due to the nature of Jeff — 1/2 wavefunction, 
Sr3lr207 bilayer material exhibits a nearly band insula- 
tor distinctly different from Sr2lr04. This originates from 
large bilayer hopping terms and the alternating rotations 
of the local octahedra. We then study the effects of elec- 
tronic interactions using a multi-orbital Hubbard model 
with Hund's coupling in Sec. |III| A first order transition 
is found within mean field theory, and the jump in the 
magnetization and the critical interaction strength de- 
pend on the SOC value. The implications of our study 
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in the context of recent experimental results are discussed 
in the section Hvl 



II. TIGHT-BINDING MODEL 
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FIG. 1. [color online] Tight binding hopping parameters: the 
red and blue dots represent two difterent Ir atom environ- 
ments due to staggered octahedra rotation. {t,ti,tz,ts,tl) 
are NN hopping between the same orbitals, while it' ,t'^) be- 
tween different orbitals. (t„,tzn,t„d,tti)/{'tid) are next NN 
hopping integrals between the same/different orbitals. 

The bilayer structure is represented by a dashed cube 
in Fig. [l] with the red and blue vertices denoting the 
two different Ir atoms in the bilayer structure due to 
IrOg staggered rotatiorPSl (the in-plane unit cell area is 
thus doubled). Taking into account the t2g d-orbital 
a = yz,xz,xy with momentum k, sublattice 7 = 
blue(B), red(R), and spin a, the tight binding Hamil- 
tonian can be written as '^k'^t i-^o i^k.i', with spinor 

V'fe.i = (<t'^<t"^<t'''''^ ^ RA^^ir where 
dl'k'f is a annihilation operator at the layer /,/' = 1,2, 
sublattice 7 and orbital a. Using this basis, the matri- 
ces Hq including nearest neighbor (NN) and next NN 
hoppings, is given by 



TjlV Tjll' 



11'^ 

BR 



W + [time-reversed], 
(1) 

with Hao being the atomic spin-orbit coupling, AL^ • S^. 
The intra-layer intra-sublattice hopping (next NN) of 
^^BB= Hrr are given by 



Hso + - 



(2) 



4i„ cos(/ca;) cos(fcj^) -I- 



where e^^ 

Atnd cos kx cos ky, and cid = 'itid sin k^ sin ky. ^xy de- 
notes the atomic potential difference between xy and one- 
dimensional xz/yz orbitals due to tetragonal distortion. 
Hgj^ is NN hopping terms between B and R sublattice 
sites, and is written as 
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where the dispersions are given by 



2{ti cos{ky) 



cos(fc,)). 



2(^1 cos(fca;) -|- t° cos{ky)), 



2i(cos(fc^) + cos{ky)), and e™* = 2t'{cos{kx) + cos(fcy)). 
For the bilayer hopping terms, it is important to notice 
that a R (B) atom on one layer lies on top of a B (R) on 
the other layer. Therefore, the bilayer hopping terms in 



i?o^ are given by 
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Itzn COS ky^ e^^ 



2t,„cosfca., and 



where e^' 

COS kx cos ky. The hopping parcinietGrs of (t, ti , i^, i', 
t'z, t„, t^n, hd, t\ ti, tnd, tfi) are shown in Fig. [l| 

Setting t as a unit, there appears to be 11 other hop- 
ping parameters, but they are not all independent of 
each other. Note that in an ideal octahedra, fixy = 0, 
t = ti = tz, tn — tzn , and = tl based on cubic symme- 
try, while t' — t'^ — without the staggered rotation of 
octahedra. Due to the IrOe rotation, the above relations 
break down and how they differ depends on the angle 
of the staggered rotation. Using Slater-Koster theorjBH 
with tdda ■ tdd^ ■ tddS = 3/2 : -1 : 1/^, and taking into 
account a distance factor of 0.9 for bilayer terms and of 
0.2 for the next NN for the exponential suppression of 
hopping parameters with distance, we found {t, ti, tz, t' , 
t'z, tn, tzn, hd, t\ ti, tnd, ti) = (-1.0, -0.94, -0.8, 0.15, 
0.36, 0.16, 0.2, 0.11, 0.27, 0.15, 0, 0) for a staggered ro- 
tation angle of ±12°. Once the ratio between tdda, tddir, 
and tdds is fixed, only the SOC A and the tetragonal 
splitting ^.xy remain independent parameters. 

The tight binding band structure is shown in Fig. [2] 
for \/t — 3 and fXxy = 0, where Jeff = 1/2 bands are in 
red, while Jeff = 3/2 bands are in blue. There are sev- 
eral important features to notice. First, J^s =1/2 and 
3/2 bands below the Fermi level are mixed along F to 
M = {■jt/2, n/2) and along F to X = (tt, 0), and the bands 
near F are Jcs = 3/2 states, not 1/2. Second, this band 
structure is similar to recently reported ARPES datcP^^ 
on Sr3lr207 and a first principle calculatiorF', but the 
bands near F h ave been misidentified as Jeff = 1/2 in 
these works.l^I^ Third, our results imply that the SOC is 
large, but not enough to fully separate Jeff =1/2 and 3/2 
bands below t he Fe rmi level, similar to theoretical stud- 
ies on Sr2lrOii212SI Qn the other hand, the unoccupied 
bands above the chemical potential are pure Jeff = 1/2. 
Thus the RIXS intensity contains more J^s =1/2 band 
contribution, because it is from a combination of unoc- 
cupied and occupied bands, while ARPES measures only 
occupied states. 

The M — X path is the reduced Brillouin zone (BZ) 
boundary due to the staggered octahedra rotation, and 
the degeneracy is protected because the in-plane poten- 
tial generated by the staggered rotation has the form of 
cos kx + cos ky which is absent along this path. The en- 
ergy difference in Jgff =1/2 bands is largest at X point, 
because of a constructive combination of t, tz, t' , and 
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FIG. 2. [color online] (a) Underlying band structure of 
Sr3lr207 along the path shown in (c). The weight of the 
Jeff =1/2 state is colored in red while the Jeff = 3/2 weight 
is in blue, and we see that the relevant bands near the Fermi 
level is mostly composed of Jeff = 1/2 wavefunction, but note 
that the band immediately below the Fermi level around the 
r point is made of Jeff = 3/2 state. The band structure is 
distinct from the single layered Sr2lr04 system shown in (b) 
mainly due to the large bilayer hoppings inherited from the 
nature of the Jeff = 1/2 wavefunction. 



t'j,. A typical splitting of 0.8i could be measured by the 
separation of the two unoccupied bands of JcS ~ l/l 
near the F point (denoted by the blue arrow in Fig. [2^). 
In Ref. [TS] the bilayer splitting was estimated by the 
separation between the two highest occupied bands at 
r. However, the energy difference between either the 1st 
and 2nd or the 1st and 3rd highest occupied band at 
the F point in Fig. [2f^a) is due to both the bilayer hop- 
pings and the hybridization between Jeff = 1/2 and 3/2 
bands. This splitting further gives an estimate of our 
unit t ~ 200 meV, which is consistent with the overall 
bandwidth when comparing our tight binding spectrum 
with the ARPES datfpl. 

Last, the most important result is the topology of the 
band structure. It is almost a band insulator where the 
Fermi level barely touches the bottom and top of unoccu- 
pied and occupied bands near M and X respectively. A 
small difference in the hopping parameters could change 
particular features such as the size of the small hole and 
electron pockets, but the separation of the two unoccu- 
pied bands from the rest via a direct gap at every fc-point, 
does not depend on the details. The direct gap between 
these two unoccupied bands occurs due to a finite t' or 
t'^, shown as the red arrow in Fig. [2|^a), indicating the 



importance of staggered octahedra rotation between the 
NN Ir atoms. For comparison, we present the underlying 
band structure of the single layer iridates in Fig. [2jb) 
that shows the large Fermi surface crossing along the M- 
X path. Now that we are equipped with the proper tight 
binding model, let us move to a magnetic ordering mech- 
anism. 



III. MEAN-FIELD THEORY 

A general interaction Hamiltonian in multi-orbital sys- 
tems is given by 

Hint " U n^a^nia^ + U' n^anip — J ^ S,jQ-Si^, 

(5) 

where the Hund's coupling J = {U — U')/2 is determined 
by intra- and inte r-orbi tal Hubbard U and U' and where 
we set U' = O.SlP^^. We treat these interactions at 
the mean-field (MF) level in the magnetic channels to 
find possible magnetic orderings. To consider the strong 
SOC, we define the order parameter in the JoS-basis. To 
do so, we first rewrite the above Hint in the Jeff-basis, 
and decouple all the terms in the magnetic channels. 
The mean field Hamiltonian is then given by Hf^f = 

~ (y + ^ ~ ¥") niil'Sil- + ^ ~ i) J2n=2,3 "1™' 

Sin. Here mi„ = (S„) = 5(c-„^(J^,.Cj„^), where <t is 
the Pauli matrix for pseudospin for n = 1,2,3 that cor- 
respond to (Jeff, Joff) = (1/2, ±1/2), (3/2, ±3/2), and 
(3/2, ±1/2), respectively. The order parameters m, „ are 
then determined self-consistently. There are 36 MF order 
parameters from three pseudospin states (n) , two sublat- 
tices (7), two layers (l), and three directions of m. 

As U increases, there is a weakly first order phase tran- 
sition from a barely band insulator into a magnetically 
ordered phase as shown in Fig. [3j mi is finite and the 
magnetic pattern is given by a G-type antiferromagnetic 
order where the moments in the blue and red atoms point 
in opposite directions. The 10.2 and m3 order parameters 
are smaller by two orders of magnitude, confirming that 
the moments are made of JcS = 1/2 electrons. The crit- 
ical interaction Uc depends on the SOC strength. The 
bigger the SOC, the smaller the Uc is, similar to the sin- 
gle layer case reported in Ref. and confirmed here, 
but opposite to three dimensional materials, SrlrOgP^. 

A finite moment affects the band structure, where the 
bottom of electron and top of hole bands are pushed away 
from the Fermi level as shown in Fig. [4j making a charge 
gap visible. For comparison, we show the band struc- 
ture after magnetic ordering for the Sr2lr04 system also, 
where the bands are much more affected by the order- 
ing than in the Sr3lr207 system. When U > Uc ^ it 
for \/t — 3, there are two competing ordering patterns 
- canted antiferromagnet (AF) and collinear AF. We 
found that when the ratio of t'^/tz is equal to the ratio of 
t' /t, the canted AF and collinear AF become degenerate. 
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FIG. 3. [color online] First order phase transition is found for 
bilayer (BL) Sr3lr207 for \/t — 2.5, 3 and 3.5. For compari- 
son, the transition for single layer (SL) Sr2lr04 is also shown 
for the same SOC strengths. We note that for the same set 
of parameters, the Sr3lr207 requires a higher Uc for the mag- 
netic moment to set in, due to its underlying band insulator. 
Note that only |mij is shown in the figure since |m2| and Img] 
are negligible. 
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FIG. 4. [color online] Band structure with magnetic ordering 
for Sr3lr207 (and Sr2lr04 in inset) using the self-consistent 
solution at U/t = 4. The charge gap from these calculation 
are Ac ~ OAt for Sr3lr207 and Ac ~ 0.6t for Sr2lr04. 

From the Slater-Koster theory shown above, t'^/tz « 2.4 
t' /t, and a G-type collinear AF state along the c-axis is 
favored by the bilayer coupling. Such a state has been 
recently confirmed in experimentJ^^''^. 

Given the small magnetic moment observed in 
Sr3lr207^, the system is close to a magnetic transition, 
where the size of the charge gap from transport data does 
not support a strongly correlated insulat or .^21 Indeed when 
the magnetic moment disappears, the band structure is 
almost a band insulator where the Fermi level barely 
touches the bottom/top of electron/hole bands. How- 
ever, the magnetic moment appears only in Jeff = 1/2 



channel, and one may ask if this magnetic insulator is 
connected to the large J7-limit of spin-orbit Mott insu- 
lator. The spin-model derived from a Jeff — 1/2 only 
model"^^ indeed displays the same magnetic order, sug- 
gesting that the two limits are adiabatically connected. 



IV. DISCUSSION AND SUMMARY 

Recently, various measurements on Sr3lr207 us- 
ing different techniqu es, suc h as ARPES^^^^, neutron 
scattering, and Rixg 9 | i4 | 32 | addition to transportP^Sl 

and optical studied, have been reported. It is important 
to ask whether the current microscopic theory is compat- 
ible with the measured quantities, and further, if it offers 
a useful starting point for future studies. 

First of all, the latest ARPES datgPI^Sl provide the dis- 
persion of the occupied bands. Using Slater-Koster hop- 
ping integrals for d-orbitals and considering the local dis- 
tortion of the octahedra as well as a distance factor that 
is exponentially decaying with distance, we obtained the 
hopping parameters up to next-nearest-neighbor. Note 
that the only independent parameters left were the SOC 
strength and the tetragonal distortion. Comparing the 
computed band structures and ARPES data, we show 
that the band below the Fermi energy is indeed mainly 
Jeff = 1/2, except near the F point where there is a large 
contribution from Jeff = 3/2. Overall, our tight-binding 
band structure fits well wit h first principle calculations^ 
and the recent ARPES datcP^l 

Furthermore, the current mean field study shows that 
the jump in the magnetic moment at the first order tran- 
sition is weaker in the bilayer iridates than the single layer 
iridates. This is due to the band structure being nearly 
insulating, advancing our understanding of the smaller 
magnetic moment in bilayer iridated^. As a consequence, 
the size of the charge gap is smaller in the bilayer sys- 
tem com pare d with the single layer one, agreeing with 
transporlPi^ and optical conduct ivitJ^^ data. 

A remaining puzzle is the giant magnon gap of over 
90 meV reported in the RIXS datfP. Our microscopic 
model contains terms that break the spin rotation sym- 
metry, a source for the magnon gap. Using a second 
order perturbation theory, this microscopic model yields 
a spin-wave spectrum with a magnon gap of approxi- 
mately 35 meV, assuming that t — 200 meV and U/t — 4. 
The gap originates from the c-axis anisotropic exchange 
terms being different from the in-plane ones. Since the 
tight binding parameters fit the ARPES data well, this 
discrepancy in the size of the magnon gap implies that 
the second order perturbation assuming a large U/t limit 
is not appropriate to estimate the magnon gap reported 
in the RIXS spectrum, and one should add higher order 
terms in the perturbation theory. In this intermediate-C/ 
range, the spin susceptibility using the random phase ap- 
proximation might be a more appropriate approach than 
the semi-classical spin wave theory applied to the spin 
model obtained in the large-?7 limit. Our microscopic 
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model is a useful starting point for future study of the 
spin susceptibility for general wavevector and frequency. 

In summary, we build a tight binding model for 
Sr3lr207 and show that the non-interacting system is al- 
most a band insulator where the two unoccupied Jeff = 
1/2 bands are separated from the rest of occupied bands 
by a direct gap at every k-point. As the interaction 
strength increases, a band insulator to magnetic insulator 
transition occurs, and a finite magnetic moment pushes 
the top/bottom of hole/electron bands further away from 
the Fermi level making the direct band gap bigger. This 
is qualitatively different from the single layer Sr2lr04, 
where the degeneracy of the bands crossing the Fermi 
level aloiig M — X is lifted turning it into a magnetic 
insulator. ■ The ground state magnetic ordering pattern 



is sensitive to the lattice structure, with the staggered 
rotation of IrOe octahedra between adjacent layers play- 
ing a crucial role in both developing a band insulator in 
the tight binding spectrum and determining the canted 
AF ordering pattern. Sr3lr207 is a magnetic insulator 
with a small momenif^J, and thus it is close to the tran- 
sition. However, the AF ordering pattern obtained from 
the spin model derived in the large J7-limit'^^ is identical 
to that obtained from the current mean-field theory, im- 
plying that this small moment insulator is likely smoothly 
connected to the Mott insulating regime as U increases, 
where the charge gap persists above the magnetic order- 
ing temperature. 
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